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FOREWORD 


This Final Report describes the efforts on Task IV during the period 22 August 1969 to 31 April 
1970. The first three tasks under this contract were concerned with the collection and assessment of 
ideas for the exploitation of the space environment. These tasks were especially concerned with. 

I Growth of high density, highly perfect crystals. 

II. Melting, purification and solidification of materials 

III. Unit separation processes. 

Work on these first three tasks is monitored by Dr. E C McKannon, Marshall Space Fhght Center, 
and is reported separately. The work reported herein is concerned with a fourth task 

IV. Field management for positioning and processmg of free suspended liquid materials. 

The contract monitor is Leroy H. Berge, Marshall Space Fhght Center. This work has as its objective 
the definition in some depth of the physical problems associated with the melting, purification, and 
solidification of free floating materials m a weightless environment. Another objective of this Task 
is to define experimental techniques and hardware which will be required in developing new 
metallurgical and ceramic processes for exploitation of the space environment. 
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ABSTRACT 


Many new processes which can exploit the weightless environment of space have been suggested as 
possibilities for making improved or unique materials. A large number of these will involve, at some 
stage, a containerless melt or transfer of molten material. Some of the physics and technology 
problems associated with these processes are discussed. The range of applicability of some new 
electromagnetic process control methods is also presented. Development of a servo device for 
measurement and control of free floating masses is described. 
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I. 


INTRODUCTION 


The prospect of the availability of manned experimental facilities in the nearly weightless 
environment of near earth orbit opens up a number of exciting possibilities for basic studies in 
metallurgical and ceramic processing. Even at this point in time it appears likely thatsome of these 
experiments can lead to basic new insights into phenomena related to materials processing methods 
and, hopefully, to future commercial exploitation of new methods for production of selected 
products of high value. As space transportation costs are reduced with development of the space 
shuttle and reuseability concepts, the range of feasible materials processes will increase rapidly as 
products of lower dollar value per pound can be considered. It goes without saying that initial 
experiments should emphasize heavily the acquisition of new basic knowledge and that for some 
time m the future only commercial products having an extremely high monetary worth per unit 
weight can be considered for pilot production studies. 
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n. 


EXAMPLES OF WEIGHTLESS MATERIALS PROCESSES 


In two symposia which NASA has held at the Marshall Space Flight Center, a large number of 
specific suggestions have been made for new types of processes or for improvements in existing 
materials preparation methods which may be applicable m the future to commercial production of 
uniquely valuable products of relatively modest weight. This report will not attempt to summarize 
these ideas for exploitation of the removal of weight forces. The reader is referred instead to the 
published proceedings of these symposia, which also mclude many other references to zero gravity 
space processing suggestions, and related work. We will here summarize only some suggested new 
processes which make use of the following categories of phenomena attending the removal of 
weight forces 

1. The elimination of buoyant separation of melt phases of different densities 

2. The elimination of thermal convection currents. 

3. The possibility for crucibleless melting of metals and ceramics. 

4. The possibility for solidification m the absence of contact with molds. 

5. Shape forming with surface tension and electromagnetic field forces. 

We can think of many processes which exploit these new possibilities either singly or in 
combination. 

A. Elimination of Gravitational Phase Separation 

Under the first category, namely the elimination of gravitational separation of melt phases of 
different densities, we can consider the following candidate products, which are a partial listing of 
ideas suggested to date - 

1 . The preparation of new alloys from immiscible melts. 

2. The preparation of mixtures of metals with metal oxides, carbides, nitrides, etc. 

3. The preparation of foams of metals or glasses. 

The occurrence of two-liquid immiscibility in the hypermonotectic alloys at temperatures above the 
monotectic introduces structural considerations. Given time enough, the two liquids will separate 
into two layers placed according to density, with the lighter layer water on top. It is quite possible, 
however, to have the two liquids existing as an emulsion with tiny droplets of one liquid suspended 
in the other liquid. Which of the two situations, the layer or the emulsion, will be obtained during a 
given heating and cooling cycle will depend upon the physical characteristics of the alloy system 
concerned, upon the conditions of formation of the second liquid, and upon the opportunity 
afforded for the segregation of the two liquids. The possibility of circumventing the buoyant 
segregation of the two liquids m the weightless environment offers the possibility of preparing 
mixtures of these melts with controlled concentration of the mixture of solid phases. 
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A well know example of an alloy system exhibiting liquid immicibility beyond the monotectic 
point at suffificiently high temperatures is the Cu-Pb alloy system (II-l). For weight percentages of 
lead above the monotectic composition point, cooling from sufficiently high temperature will cause 
breakdown into two liquid layers. The Cu-rich layer forms the top layer and the Pb-iich layer the 
• bottom layer. Some interesting systems are systems such aCs Cu-W which are virtually insoluble in 
[ each other both in the liquid and solid state. One again gets a two layer system and the bottom 
layer W solidifies first at the freezing point of W, the top layer freezing at the freezing point of Cu. 
Of interest would be the preparation of an emulsion of W and Cu which would be readily 
. accomplished in the weightless state and examination of the properties of the mixture produced 
, upon solidification. 

The preparation of many mixtures of metals with metal oxides, carbides, nitrides, etc., involves 
components with widely different densities. The solubility of oxides, nitrides, etc. in metals is often 
measured in parts per million. Vigorous stirring will cause mixing but within minutes or even 
seconds there is separation due to density differences. In a weightless situation, there would be no 
tendency to segregate due to density differences. Electromagnetic positioning devices discussed in 
this report would also aid in stirring the components to achieve the dispersion qualities desired. The 
preparation of dispersions of many substances not contemplated on the Earth could be considered 
here, such as unique cermet-like \materials for structural and electronic applications. 

Some examples of these mixtures which have been suggested are rare earth oxide dispersions m 
superalloys (II-2), oxides of fissionable materials in metallic reactor fuel elements, and dispersions 
of small particles of glass in a matrix having a different index of refraction to form phototropic 
materials or Christiansen filters. However, with the possibility of producing many different mixtures 
as discussed above, there is a unique opportunity for developing totally new materials not 
contemplated by present technology. 

The possibility of producing foams of metals or glasses in the weightless environment offers the 
possibility of producing very light materials with good service properties. One of the major 
problems m foam production is the problem of drainage. In the weightless operation, the problem 
of foam drainage will be eliminated, even though other problems will remain, such as the collapse of 
bubbles due to surface tension forces or drainage in thin films through suction by Plateau’s border. 

Some suggestions for application of foam materials are crushable structures, battery electrodes, the 
production of very light materials with high temperature capability, high strength properties, and 
immitation of natural products like wood (II-3). Wood has good properties under compression in 
that it is composed of cells. When the yield point is reached the cells are crushed upon further 
application of load. This reduces the buckling of these materials. The design of materials with cell 
like structure would be a very useful imitation of nature. 
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B. 


Elimination of Thermal Convection 


The effect of thermal convection in forming dislocations m crystals grown from melts has been 
described by Utech (II-4) and Grodzka (II-5). Hamalainen (II-6) has also described frozen Benard 
convection cells in thm alkali halide crystals. This work demonstrates the importance of thermal 
convection upon the crystal structure obtained. The effect of thermal convection dunng glass 
formation on the homogeneity of glasses has been discussed by Deeg (II-7). It seems likely that, in 
many other cases, the effects of thermal convection may have as yet unrecognized effects. 

C. Crucibleless Melting 

Limited studies of the potential for crucibleless melting of metals have been performed in terrestrial 
experiments where small spherules of metals have been levitated and melted by use of radio 
frequency induction fields (II-8 through 11-22). The objectives of this work have generally fallen 
under one of the following categories 

1. Melting and reaction of refractory or reactive metals, some of which are 
contaminated by any contact with crucible. 

2 Formation of an alloy completely free from segregations. 

3. Observation of subcooling below the normal freezing point, which can be achieved 
routinely with some metals. 

The method has not found commercial application for several reasons. First, levitation is not 
practicable for poor conductors. Secondly, the required levitation power is often so high as to 
preclude separate control of specimen temperature. In other words, a reduction in RF heating of 
the specimen will cause loss of levitation or pouring from the bottom of the specimen. Instabilities 
are sometimes encountered, and, in any event, the mass of material which can be prepared is 
severely limited. Since levitation is the most obvious property of the weightless space environment, 
we can for the first time consider a whole range of new processes of this type and their potential 
application to preparation of large quantities of material. 

Removal of the requirements for levitation of a crucibleless melt means that we can have complete 
latitude in heatmg of the specimen and its temperature control during solidification. Materials can 
be prepared which are nonconductors, and many suggestions have been made for new glass 
processes which should be possible in the space environment. 

Besides the possibility for reacting and melting materials without crucible contact, a corollary is the 
ability to provide super heating to reactive and high melting metals for which skull melting 
techmques must currently be used. This opens up the possibility for precision casting of these 
materials which cannot at present be accomplished with the negligible superheat available in skull 
meltmg. 
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D. The Possibility for Solidification in the Absence of Contact with Molds 

In terrestrial levitation experiments, extreme degrees of subcooling before solidification are often 
observed (11-23). Walker has achieved subcooling for a number of materials in glass-lined crucibles. 
In some cases, the degree of sub cooling is believed to approach that at which homogenous 
nucleation occurs. The possibility exists that some materials may form new phases if solidified in 
the absence of heterogenous nucleation caused by crucible or mold contact. For example, a new 
phase of gallium, denoted as Gallium 3, has recently been prepared in minute quantities at a 
temperature of -30° centigrade (11-24). 

The possibility for forming glasses from materials which normally crystallize through heterogenous 
nucleation has been discussed in papers by Olsen (11-25) and Happe (11-26). Since the nature of such 
products cannot be anticipated at present, it is clear that containerless solidification experiments 
must initially be proposed with emphasis upon obtaining basic information rather than with any 
attempt to anticipate commercial applications. 

It has been suggested by Witt (11-27) that the difficulty encountered in pulling thm filaments of 
single crystal material from melts in studies by Monsanto Research would be greatly relieved if it 
were possible to suspend a crucibleless melt of boron, for example, and apply forces to the melt 
which would oppose forces due to withdrawing a boron fiber at a rate allowing single crystal 
formation. 

The possibilities for formation of near perfect spheres by solidification of melts under the sole 
action of their own surface tension has also been suggested. 

E. Shape Formation Through Surface Tension and Inertial and Electromagnetic Fields 

Deeg (11-28) has suggested possibilities of fire polished glass surfaces for spherical or spheroidal 
surfaces which could be formed by centrifugal action in a freely suspended rotating melt. Rotations 
to produce an oblate spheroidal form of given eccentricity can be easily imparted to a floating mass 
by means of the orthogonal coil sets discussed in a following section. 

Alternatively or in conjunction with inertial forces produced by rotation, electrostatic fields can be 
applied to perturb the spheroidal shape. If the electric field at the surface of a conductor is 
represented by E, the local curvature imparted to the surface will be given by 



where r j and T 2 are the radii of curvature in two orthogonal directions, p is the internal pressure, a 
is the surface tension, and tc is the dielectric constant. Where only slight departures from spheroidal 
form are required, design of suitable electrodes to furnish fields to produce a desired shape may be 
feasible. ■ 
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F. Positioning and Handling of Containerless Melts - Theoretical and Experimental Studies 

A large number of the suggested new processes in weightless processing of metals and ceramics 
involve at some stage a floating molten mass out of contact with crucibles or molds. A floating mass 
of this type, held together solely by its own surface tension force and free to oscillate in shape, 
rotate and support internal fluid currents which might be driven by these or other applied forces, 
represents a phenomenon with which we have had little experience. Such phenomena are briefly 
encountered in shot tower processes and in free falling liquid streams, but here the effects of air 
resistance and, usually, lack of time to establish equilibrium exist. 

The following section outlines preliminary results' on investigation of some physical problems 
associated with the melting, handling, and solidification of such freely floating masses. The study of 
the basic physics involved in such processes leads quite naturally to a definition of the experimental 
facilities which will be required to handle a range of suggested processes. This section also describes 
some preliminary development of some new physical hardware concepts which will be required in 
such facilities. 
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HI. OSCILLATIONS OF A FREE LIQUID SPHEROID 


The mathematical physics of small shape oscillations in a liquid mass under the sole action of 
surface tension and viscous and inertial forces was treated by Lord Rayleigh (III-l), and more 
recently by Chandrasekhar (III-2).The extension of this work to include large oscillations and large 
viscosity have been made by Benedikt (III-3) and Reid (III-4). For reasonably small oscillations. 
Figure III-l gives the oscillation frequency and time decay constant for shape oscillations in which 
the shape is alternately a prolate and oblate spheroid. This is the simplest mode of oscillation. More 
complex shape oscillations can also occur; for a given material the oscillation frequency will be 
greater and the decay time shorter than for the mode for which Figure III-l applies. The 
mathematical expressions for the frequency and decay tune for a given oscillation mode are 
available from the classical literature (III-5). 

Figure III-2 illustrates the first three hydrodynamic modes. We can describe these modes most easily 
by expressing the distance from the center of mass to the liquid surface r as a function of colatitude 
0 in the way shown at the bottom of the figure. For the zero mode, the radius vector equals a 
constant a Q multiplied by the 0-order Lengendre polynomial which is simply unity. For this mode, 
the mass is a perfect sphere with no distortion. For the second hydrodynamic mode, which 
represents a shape oscillation between prolate and oblate spheroid, we must add to the original 
radius a term ^ 2^2 ( cos cos 0 } 2 *' The cen ^ er drawing in Figure III-2 illustrates this shape at a 
given instant of time. The shape oscillates sinusoidally between oblate and prolate. This oscillation 
mode can be either excited or damped by application of the magnetostriction forces through 
excitation of position control coils. 
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Figure III-2. Hydrodynamic Modes 


A. Mass Limitation 

Mass limitations arising from excitation of oscillations are described in Reference III-6. Liquid metal 
masses at least as large as tens of kilograms under control accelerations of 10'^g can be considered 
without worrying about rupture of the mass due to the necessity of position control in practical 
cases. 

As discussed in the following section, studies have been made of methods for position control of a 
free floating spheroid of liquid. The method which will be described in detail involves the use of 
electromagnetic fields provided by three orthogonal sets of coils excited by high frequency currents. 
Regardless of the method used for position control for such free floating liquid masses, it can be 
anticipated that excitation of shape oscillations and internal fluid currents will occur. 

If the position control coils along a given axis are excited asymmetrically, e.g., if we excite only 
one, then we will obtain a net translational force which excites both the first and third mode. The 
third mode resembles an eggplant and can be treated similarly to the second mode. The first mode, 
shown as P j in Figure III-2 represents pure translation of a spherical shape and must be treated 
differently from the other modes. The reason for this is that the frequency corresponding to the 
first mode is zero and, for this particular term in the sum, we should replace the sinusoidal time 
term by a term lmear in time. 

B. Free Oscillations 

The free oscillations of a viscous liquid globe under the action of surface tension forces alone are of 
interest in determining the natural frequencies of the system where resonances can occur. The basic 
problem we consider here is a fluid sphere of incompressible, viscous fluid under the action of 
surface tension alone and the small oscillations of the fluid about its spherical form. Lamb (III-5) 
had discussed this problem in detail for inviscous fluids. Smce we are considering small departures 
from an equilibrium spherical state of an incompressible fluid, we can describe the deformed surface 
of the configuration by the equation 
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r = R [ 1 + € Y (0, $) 3 
for the deformations of the liquid surface. 

The function e is a function of time and here the dependence of the velocities and other quantities 
on time enters explicitly. It is supposed that 

-O' t 

e = e e 
o 

where e 0 is a constant and a is composed of a real part expressing the damping of oscillation due to 
viscous forces and an imaginary part expressing the oscillatory nature of the deformations about the 
spherical form 

-O' t -t/T -i27rft 
e = e e 

where r is the damping time and f is the frequency of the oscillation. 

The angular frequencies of oscillation are given by 

(2irf) 2 = o i ;c 2 = A (JL - 1) (JL + 2) Tj/p R 3 

where 

a ^ . Q = the angular oscillation frequency for inviscous liquids 
J t ~ the degree of the surface harmonic Y™ (0, <i>) 

R = the radius of the undeformed sphere 

p — density of the fluid 

T i = surface tension/unit length 

Y P^(cosfl) jg^cosm^ + h^sinm^J 
where the P™ are the associated Lengendre polynommals. 

For small oscillations of viscous fluids about their spherical form, the problem is to determine the 
0 £ * 8 or the frequencies of oscillation for a coefficient of kinetic viscosity v. This involves solving 

the linearized equations of motion' 
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5 u 
a t 



- v V x (V x u) 


V- u = 0 


where u = velocity, P = pressure, p = density, and t is time, subject to the following prescribed 
boundary conditions: 

1. At the deformed surface of the drop, the radial component of the velocity must be 
equal to the velocity of the surface itself. 

2. The tangential stresses vanish at the surface. 

3. The normal component of the viscous stress is -P^ = T j (1/R j + l/R^). Here R j and 
R 2 are the principal radii of curvature taken positive when directed inward. 

Reid (III-4) has discussed this problem and he obtains the characteristic equation for the 
frequencies of oscillation: 
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where C) ^(q) is written in terms of spherical Bessel functions 
Here we have 
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We can solve the characteristic equation and get the a 4 ^,’s. The following results are obtained: 
As q-* 00 , the limiting case of small viscosity 

0 *,„“<* + 1) iVR 2 ± ‘er t ;0 

which is Lamb’s result. 
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2. The characteristic equation admits solutions of two types, depending on a value of 
a called the critical value a c . 

a. If oP > a < p then damped oscillations will occur. 

b. If oP < a ( p ) two aperiodic modes of decay appear. 

Figure III-3, from Chandrasekhar’s book, illustrates the aperiodic decay constants in terms of the 
undamped angular frequency a % 0 and the physical properties of the hquid spheroid. 

C. Excitation or Damping by Positioning Device 

Excitation of the position control coils can be used either to initiate or to damp translational 
motions of the floating mass and shape oscillations as well as any rotations which may occur. The 
position control and shape oscillation problem is most easily treated by considering individual 
hydrodynamic modes of a liquid sphere under the action of surface tension forces and the 
electromagnetic forces due to internal eddy currents. 



° i, o v 

Figure III-3. Damping Constants for the Aperiodic Modes 
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D. Internal Fluid Currents 

In addition to generation or damping of surface oscillations and translational motions, we can 
stimulate circulating currents within the floating mass by magnetostrictive forces. Figure III-4 shows 
a simple current pattern which can be stimulated by exciting a single control axis. These fluid 
currents can be oriented in various directions by exciting various sets of coils. These currents will be 
useful for stirring of alloys or dispersions and will also affect bubble motion in an important way. 

At the higher frequencies, the pressure gradient resulting from the magnetostrictive forces which 
drive the steady state fluid currents will be concentrated in a thin surface layer centered about the 
equatorial plane with respect to the excitation axis. After the coils are excited, an equilibrium will 
be quickly reached in which this driving pressure gradient is balanced by viscous forces. At lower 
frequencies where the skin depth becomes an appreciable fraction of the sphere radius, the flow 
patterns will be altered somewhat. The quasiequilibrium shape of this spheroid will be influenced 
not only by the magnetostrictive forces which drive the fluid currents, but also by the dynamics of 
the circulating dense fluid Since magnetostrictive forces can be applied in an arbitrary manner and 
translational acceleration of the mass can be suppressed by simultaneous excitation of opposite 
members in coil pairs, we see that considerable freedom exists for electromagnetic shaping of the 
floating mass in various symmetric or asymmetric configurations. 

E. Rotations 

If we consider an orthogonal pair of excitation coils to be excited with a phase quandrature so as to 
produce a rotating magnetic field, we can calculate that the spheroid will be caused to rotate in the 
manner of the rotor of an induction motor. By proper choice of excitation for various pairs of the 
positioning coils, angular velocity can be imparted to the floating spheroid m any desired direction 
and with any speed consistent with configurational stability of the mass This will be limited by the 
stability furnished by the action of surface tension. Oscillational instability can occur for higher 
rotational speeds, and damping due to viscosity will play a key role m prevention of buildup of 
catastrophic oscillations. Conversely, properly phased excitation of a given pair of coils can be used 
to remove an angular velocity orthogonal to both coil axes if such rotation is not desired. 

F. Sound Wave Generation 

The high frequency variation of magnetostriction forces can be used to generate sound waves within 
the floating mass if we desire to go to high level excitation. This might be useful, for example, to 
prevent agglomeration of oxide dispersion particles within the mass. 
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Figure III-4 Fluid Currents Due to Magnetostriction 
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IV. THEORY OF ELECTROMAGNETIC POSITION CONTROL AND BODY DISTRIBUTED 
CURRENTS AND FORCES 

The present section describes theoretical developments pursued in order to examine a few idealized 
basic problems which have immediately applicable results and which can be logically extended to 
obtain solutions to other problems. The objectives of these theoretical development are- (1) to 
understand the electric currents induced in material bodies of selected composition when these 
bodies are immersed in electromagnetic fields, (2) to understand the forces exerted on freely 
suspended material bodies by electromagnetic fields, and (3) to determine the interaction of 
material bodies with the electromagnetic field sources through mutual inductance and capacitance 
so as to describe how the positioning and control devices must function. 

A. RF Position Control of Levitated Mass 

For three-axis position control of the levitated specimen in zero-g environment, three pairs of 
centering coils can be used, each pair with an axis orthogonal to the other two. The theory of the 
forces delivered to the levitated mass resultmg from excitation of any one of the coils has been 
studied by Okress (IV-1). The theory has been generalized for multiple turn coils by Fromm and 
Jehn (IV-2). Since relatively small forces are required for centering, it is expected that quite low RF 
powers will suffice. 

Figure IV-1 (from Fromm and Jehn) shows a cross section of a levitation coil used m present one-g 
experiments. In a zero-g experiment, one can utilize relatively weak currents in a few turns and 
provide stability by the addition of a separately powered coil lying above the sphere, and similarly 
other sets of coils in a position, to impart forces in all three orthogonal directions. Excitation of 
any one of the coils would cause a repulsion away from it so that a total of six coils in three 
orthogonal sets would do thej'ob. 

Figure IV-1 shows the way m which the force of repulsion depends on the dimensionless parameter 
X, the ratio of sphere radius to skin depth. The restoring forces on the levitated object will be 
concentrated within a “skin depth” given by the formula in the figure. Here, f is the frequency, <y is 
the conductivity, and n is the permeability. Larger forces are available, at a given power, by the use 
of higher frequency or lower resistivity materials. To some extent, the distribution of the body 
forces can be adjusted by adjustment of the frequency, m the sense that the skm depth may be 
made to penetrate through the body of the specimen by going to a sufficiently low frequency. This 
would allow imparting of restoring forces with minimum distortion of the surface shape of the 
specimen. If separate facilities are provided for resistance oven or RF heating, separate control could 
be achieved for melting and stirring, the latter achieved by use of the high intensity RF fields. 

B. Vector Potential Induced in Conducting Sphere by Uniform Oscillating Magnetic Field 

The first problem considered is that of a conducting sphere in a uniform time harmonic magnetic 
field. The solution of this problem is applicable to the study of melting conducting materials by 
induction, fluid currents produced by electromagnetic forces, fluid oscillations due to 
electromagnetic forces, studies of bubble motion driven by fluid currents, the generation of acoustic 
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waves, and many other topics of interest. This solution is also the first step in understanding 
positioning forces exerted on the sphere by nonuniform fields, since it can be used to construct a 
first-order solution for force exerted by fields having gradients, as will be shown in following 
sections. 

The problem considered is a ^sphere of uniform conductivity a immersed in a time harmonic 
uniform magnetic field B^= B Q k cos cot, where B 0 is the amplitude of the applied field, v - cof2ir 
is the frequency of the applied field, and k Q is the unit vector in the z direction. Mathematically the 
problem can be formulated as the following boundary value problem: 

V z A g = 0 outside the sphere 

V 2 iA r p(7d Aj/S t = 0 inside the sphere \ 

and boundary conditions: 

As r^oo , Aq-* A 

Asr-* o, Aj remains finite 

At r = a, n * [( V x Ap — (V x A e )] = 0 

n x [ 1/Mj v x Aj - l/p 0 V x A e ] = 0. 

Here A is the vector potential and the magnetic induction B is just V x A. 

Thus 


B^ = V x A^ is the applied magnetic field 
B e = V x A e is the total external magnetic field 

Bj = V x A^ is the magnetic field within the sphere, and solving the boundary value 
posed above for A e and A_j gives B e and Bj. 

The above problem has been treated by Smythe (11-20) and others. We have extended the treatment 
of Symthe to include the field inside the sphere, the instantaneous forces acting at all points on and 
within the sphere, and the time averaged forces acting at all points on and within the sphere. Rather 
than present the detailed solution to this problem which can be found in the above reference, we 
shall outline the solution and discuss the results and our extension in detail. An alternate solution m 
terms of the co - variant form of field equations is given in Appendix G. 

To simplify the mathematical solution of the problem, it is convenient to replace B ^ = B Q cos cot k 
by = B 0 ei w *k, where = cos cot + j sin cot, where j = 
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We can then write 


A - j a) M cr 0« 

We must remember that in the solutions we must multiply by e^* and take the real part as the 
actual solution. 


The vector potential of can be written in spherical polar coordinates as A A = Vi B Q r sin 0 <&, 
where & is the unit vector in the $ direction. Since the exciting vector potential has only a $ 
component, we expect the currents induced in the sphere to have only a <t> component, and hence 
the vector potential Aj to have only a $> component. Since the total external vector potential is due 
to the applied field and the induced A currents, we then have the vector potential external to the 
sphere; must also have only a & component. The results obtained by S my the (11-20) after 
writing V ^ A=0 and V 2 A-jco/ru A = 0 in spherical polar coordinates and solving the boundary 
value problem by separation of variables are: 


e 


= [irB rsin0+DB sin0/r ]$ 


A. = h.B 


2JD C 
i 2 o r 


_ i 

2 


Ig/ 2 {(jp) 2 r}sm0& 


where 


(2 ft + M c ) V I_ 1/2 M - 

(1 + v 2 ) + 2 m 

h/2 {V 

(M -K 0 )v I. 1/2 {v} + , 

2 

M o a + v ) - m 

h/2 < V > 


and 


C = 


3/x v a 


3/2 


r 2 

0* “ f* Q ) V K 1 /2 { V } + (1 + V ) - fi I x y 2 { v} 


Here we have I n as the symbol for the Bessel function of order n,p = ofico and v = (jp)^ a. The skin 
depth 5 is a basic parameter of the solution and 


6 = 




2 

CT JU Cl) 
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If we consider the conducting sphere as not permeable so that p = p Q , then C further simplifies to 
become 


C = 


3 a 


3/2 


Vl l/2 {V} 


Explicit Analytical Expressions for Body Distributed Currents and Forces 

Eddy Currents 

An expression for the magnetic potential A and the current density j is available from the work of 
Smythe (11-20) as outlined above. 


Figures IV -2 and IV-3 show the results of calculations of the radial component of the magnetic field 
inside a conducting sphere placed in a previously uniform oscillating magnetic field. The ordmants 
give the radial component of the magnetic field inside this sphere (B^) normalized to the 
unperturbed value B of the external field. The parameters labeling the various curves are the ratio of 
sphere radius R 2 to skin depth 8 . As noted above, the skin depth is determined by the angular 
frequency co, the sphere conductivity a, and the permeability p. High frequencies and good 
conductivity correspond to small depths of penetration 8 of the electromagnetic field into the 
sphere and, correspondingly, to large values of the parameter R 2 /S. We see that for large skin depth 
(R 2/5 = 0 . 1 ), the field penetrates far inside the sphere and is, in fact, nearly uniform throughout the 
sphere. For high frequencies and good conductivity (R 2/8 =13), the field dies out rapidly inside the 
surface of the sphere and, in fact, reverses sign as we proceed inward. 


Figures IV-4 and IV-5 show similar calculations for the 8 -component of the magnetic field 
computed from curl Aj. 


The force density is obtained in terms of the curl of this expression and the current density for the 
case of unit permeability. The current density is given by 



V 2 ^ 

r* [cosh v* - — smh v* ] 

7T v 1 v 


where we have written the Bessel functions m terms of hyperbolic functions with the notation 
B Q e^ = applied field 


= OOJfl 

a - conductivity 

p = permeability 

- 

a = sphere radius 


= -yfipo 
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RADIAL D I STANCE ( CENTIMETERS) 


Figure IV-2. Radial Component of Magnetic Field per Unit Magnetic Induction in Aluminum Sphere 
of Radius R2 = 0.4125 Inches at Colatitude 6 = 45°, Showing Penetration of Field into 
Sphere as a Function of Skin Depth 
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Figure IV-3. 


® R 2 /S = 6 
© R 2 / 8 » 8 



RADIAL DISTANCE IN CENTIMETERS 


Radial Component of Magnetic Field per Unit Magnetic Induction in Aluminum 
Sphere of Radius R 2 - 0.4125 Inches at Colatitude 6 = 45°, Showing Penetration of 
Field into Sphere as a Function of Skin Depth 
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RADIAL DISTANCE ( CENTIMETERS) 


Figure IV-4. 6 Component of Magnetic Field per Unit Magnetic Induction in Aluminum Sphere 

of Radius R 2 = 0.4125 Inches at Colatitude 6 = 45°, Showing Penetration of Field 
into Sphere as a Function of Skin Depth 
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If we use the notation 



we have 

smh v = sinh if cos if + 1 cosh if sin if 
cosh v 1 = cosh if' cos if' + i sinh if' sin if' 

It is necessary to obtain the real part of both the current density and the magnetic field before 
multiplying to obtain the force density. For this it is useful to use the notation. 

t t r 

a = sinh ^ cos 4 oi = sinh 4 cos 4 7 = cosh 4 cos 4 * 

, r t r 

/? = cosh 4 sin 4 j8 = cosh 4 sin 4 17 = sinh 4 sin 4 

in terms of which the current density becomes 


1 


0 


3 B a sin 0 
0 


2 2 
2/i r 6 (a + j 3 ) 


(Y- iX) 


where we have utilized the abbreviations 
X = - 01 F + /JE 
Y = aE+(3F 

and 


E = y-T)~ 


6 

r 


a 

r 



Figure IV-6 shows results of calculations for the current density inside a conducting sphere for 
several values of skin depth. These have^ been calculated by taking the time derivative of the 
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magnetic vector potential. Since the magnetic vector potential has only an azimuthal component 
A^, these eddy currents flow in concentric circles whose centers lie along a line passing through the 
center of the sphere in the direction of the impressed magnetic field. 


We can then write the force density as 

f = R (j) x R (curl A) where R ( ) denotes “real part of . . .” 


or 



2 2 2 _ 
B a sin 0 
o 


2 , 2 A 2 

jur (a + p ) 
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d X 
5 r 


f = 

*0 2 


B 2 a 2 sin 0 cos 0 
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2 2 

+ 0 ) 
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These expressions may be useful in considering the resultant circulating fluid currents within a 
molten sphere. We note immediately that the force field is nonconservative, as might be expected 
for time dependent electromagnetic forces. 

Figure IV-7 shows results of calculations for the radial component of the body distributed force 
density for several values of skin depth We see, for example, that for a skin depth of 1/10 the 
sphere radius (curve 6), the magnetostrictive force dies off extremely rapidly as we proceed into the 
sphere. Although small radial outward forces are encountered for some regions within the sphere, 
the pressure is expected to be positive throughout the intenor of the sphere so that we need not be 
concerned with cavitation. 

Figure IV-8 shows results of similar calculations for the 6 component of the Lorentz force. 
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NOTE: POSITIVE DIRECTION TAKEN TOWARD SPHERE CENTER SO THAT 



Figure IV-7. Radial Component Force/Volume per Magnetic Induction Squared for Aluminum 
Sphere R2 = 0.4125 Inches at 6 = 45° as a Function of Skin Depth 
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Figure IV-8. 6 Component Force/Volume per Unit Magnetic Induction Squared for Aluminum 

Sphere R 2 = 0.4125 Inches at 6 = 45° as a Function of Skin Depth 
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FORCE MEASUREMENTS AND TECHNIQUES 


Force-of-repulsion measurements were made to support calculations on the fundamental physical 
phenomena and the electronic circuitry. It is important to make actual measurements m addition to 
calculations because of nonlinearities, which will be discussed later Since the electromagnetic forces 
depend upon the square of the field amplitude, the forces due to the separate positioning coils will 
not add linearly, but will depend, in a complex way, on the relative phase of excitation of the coil 
pairs, as discussed in Appendix C. These force measurements fall into two categories (1) large 
forces obtained in a general survey, and (2) small forces of the order of dynes, which are closer to 
the positioning forces actually anticipated in zero-gravity materials processing experiments. 

A. Measurement of Large Forces 

A simple, magnetically damped beam balance was modified to support any one of several spheres 
above a coil of wire, as shown schematically in Figure V-l . The coil was raised and lowered to vary 
the spacing between the coil and sphere. The sphere was relatively motionless as the balance was 
returned to a balanced condition for each position of the coil. The circuit employed to pass current 
through the coil is also shown in Figure V-l. The 10-kHz signal from the signal generator was 
amplified by a power amplifier and applied to the coil, the variable capacitor, and to the one-ohm 
resistor connected in senes The capacitor was set to tune the combination to the applied frequency 
and the resulting current was observed by reading the voltage generated across the resistor. With this 
apparatus, forces were measured for several different coil-to-sphere distances (center-to-center) for 
each of several spheres with a limiting accuracy of approximately ± 30 dynes. Materials employed 
for the spheres were copper, aluminum, Inconel, Dowmetal, and tungsten-carbide steel. Results for 
copper and aluminum at a coil current level of 1 .4 rms ampere are given m Figure V-2 with the 
forces normalized to one ampere turn. 

The same apparatus was used to obtain the data given in Figures V-3 and V-4, which are similar to 
those m the preceding figure, except that a different coil was used and data for three different 
frequencies are shown, viz. 10 kHz, 1 kHz and 100 Hz The currents used were about 0.2 rms 
ampere for 10 kHz and about 0.7 rms ampere for 1 Hz and 100 Hz. The broken lines in Figures V-3 
and V-4 are the results of calculations based upon the assumption that currents within the sphere 
are excited by a uniform oscillating magnetic field. There is some inconsistency in computing the 
net translational force for the case in which the applied field is uniform. This partially explams the 
discrepancy shown in the figures. If the field were uniform, no net force would exist. The procedure 
used m the calculation was to compute the currents m the sphere for the case of the uniform field 
and to assume that these currents were acted upon by the actual field near the coil so that a net 
force was obtained. 

B. Measurement of Small Forces 

Figure V-5 gives results of calculations over a wider range of skm depths and sphere-to-coil 
distances. We see that the force is maximum for a sphere-to-coil distance approximately 30 percent 
of the coil radius These calculations assume that the magnetic dipole moment induced in the sphere 
is the same as would be induced by a uniform field. In actuality, the dipole moment might be 
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Figure V-2. Results for Cu and A1 at a Coil Current Level of 1 4 rms Ampere 
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significantly different in positions close to the coil where the field gradient becomes severe Thus, 
we can expect that the actual position of this maximum will be shifted. 

A piece of plate glass 0.5-inch thick, 2-feet long, and 1-foot wide was equipped with three 
micrometer screw heads which permitted controlled tilting of the glass in one direction in 
increments of 3.9 x 10~ 5 radian for each 0.001 -inch division. The vernier micrometer heads allowed 
tilt measurements to one-tenth of this value or 3.9 x 1 0'^ radian. A transparent box approximately 
45 inches tall was placed on the glass plate and the same aluminum sphere as used previously was 
suspended from the top of the box by a fine wire, see Figure V-6 The same coil as used previously, 
(see Figure V-3 and V-4) was placed near the sphere to obtain the data given in Figure V-7 at ] 0 
kHz. The procedure used with the apparatus was to observe'the position of the sphere on the reticle 
of a telescope, tilt the table a small amount, and apply current sufficient to return the sphere to its 
original position. The amount of tilt and the mass of the sphere indicated the force required and the 
current in the coil was read as descnbed above. The accuracy of the force measurements with this 
technique is approximately ± 0.5 dyne, and most of this uncertainty is provided by the limit to 
which the position of the sphere may be read. Swinging of the sphere, which hampered observation 
of the sphere’s position, was eliminated by attaching a lightweight probe to the bottom of the 
sphere in line with the wire supporting it. The probe extended about 1 /8 inch into some 2,000 
centistoke silicone oil. This produced a degree of damping somewhat less than critical. By taking 
vertical “slices” of the data of Figure V-7, the information shown in Figure V-8 was obtained. 

This tilt-table apparatus was also used to measure small forces applied to spheres as they roll on an 
optically flat piece of glass. This techmque was not used extensively because the sensitivity was 
usually no greater than ± 0.5 dyne and only under ideal conditions was it as low as ± 0.05 dyne. 
The advantage of such a technique is to permit the sphere to range over a wider area, unrestrained 
by any supporting wire. A disadvantage is that a nearly perfect sphere and dust free surfaces are 
required. 

C. Measurement of Stability of Sphere Position Within Positioning Coils When Servo Loop is 
Inactivated 

Since the force between a given positioning coil and the floating conductmg mass is repulsive, the 
force fields of three orthogonal pairs of positioning coils will create a potential well centered about 
the common intersection of their axes when all coils are excited A two-dimensional analog of this 
potential well has been examined by exciting two orthogonal coil sets whose axes were horizontal. 
The conducting aluminum sphere was suspended on a long pendulum The arrangement was such 
that, with the coils deactivated, the ball was positioned at the intersection of the two coil-pair axes. 
With the coils activated, oscillations of the sphere were initiated by physically touching the ball. 
With the two members of each coil parr connected in phase or 180 degrees out of phase, a potential 
well was observed having a pincushion shape. This means that the oscillations were of lower 
amplitude when they were parallel to either coil-pair axes than when the plane of the oscillations 
made an ange of 45 degrees with both coil-pair axes. The period of the pendulum of course 
depended upon both the plane of oscillation and the amplitude and ranged between 1.5 to 2.0 
seconds. The natural frequency of the pendulum (~ 10 feet long) was 3.5 seconds. With the 
members of coil pairs connected in phase, very httle restoring force was apparent for ball positions 
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Figure V-6. Tilt Table Assembly for Measurement of Small Forces 





near the well center, and the nature of the oscillations were to some extent comparable to the case 
where the potential well would be a box with rather stiff reflecting boundaries. With the coils 
connected out of phase, the restoring forces had a “softer” nature, thus demonstrating the 
nonlinearity of force addition. It is well known that the axial field near the center of the Helmholtz 
coil pair vanes approximately with the cube of the displacement from the center position. This 
would give a force varying as the 6th power of displacement from the center and would be 
qualitatively similar to the force variation observed with the coils connected in phase. 

Measurement of small forces by a number of alternate methods was explored analytically and, to 
some extent, experimentally. A summary of this work is given m Appendix A 
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VI. POSITION SENSING AND CONTROL SERVO 


A. Introduction 

The present section describes development work on a servo device for position sensing and control 
of a free floating conducting mass. It was found that the eddy currents induced in the conducting 
sphere cause significant changes in the inductance of the coils used to generate the electromagnetic 
field. A circuit was developed to sense the change in this inductance as the ball moves relative to the 
coil. The ball position signals obtained in this way are then processed, including signal 
differentiation to detect ball velocity, and then used to control the coil excitations in a servo loop. 
Proper mixing of the ball position and velocity errors, and adjustment of gain and frequency 
characteristics in the servo loop allow for automatic restoring of the ball to a central position 
equidistant from all coils and removal of any residual velocity relative to the coils. Although, as 
noted earlier, the floating ball can be trapped in the potential well created by the coils when they 
are excited, there is essentially no damping of position oscillations. In order to provide for quiescent 
conditions in which drift of the floating mass can be corrected, damping in an automatic manner is 
considered essential. For some applications, such as the casting of perfect spheres, accurate removal 
of drift velocity will be essential so that position control forces may be removed during 
solidification. 

B. Description of Servo Circuits 

The control system hardware for position and damping is illustrated in Figures VI-1 and VI-2. 
Figure VI-3 is a block diagram showing details of the electronics employed in the loop (single axis) 
and the signal flow paths. The operation of the loop is as follows. 

The coils shown at the left of Figure VI-3 are accompanied by shunt resistors which provide a signal 
proportional to coil current. The excitation voltage and the current signals are amplified 
individually and drive zero cross-over circuits in the logic block. The output of the logic circuit is a 
series of pulses which drive the switches, one coil switching (+) reference and one coil switching (-) 
reference into two summers. The duty cycle of the switches is proportional to the phase angle 
between the current in the coil and the excitation voltage. This phase angle is related to the position 
of the ball relative to the coil, as shown in Figure VI-4. Therefore, if the ball is centered between 
the coils, the coils have identical phase shifts and the (+) and (-) reference switches have equal duty 
cycles. The resultant average seen by the error summers is zero. Also, the multiplier outputs are 
equal and both coils have identical excitation voltages from the power amplifiers. The conditions up 
to now indicate a quiescent state with no disturbances to the ball. 

Assume the ball occupies a position nearer to the coil shown in the upper channel of Figure VI-3. 
From Figure VI-4, we find this coil’s current to be shifted up in phase, and the lower coil, which is 
further from the ball, is shifted down in phase. Under these conditions, the (+) reference switch 
duty cycle decreases and the (-) reference switch cycle increases. The error in the upper channel 
becomes (+) at the summer output (the summer inverts), and the modulator increases the signal to 
the power amplifier and the upper coil. This acts to repel the ball back toward the center position. 
Even though the lower summer’s output goes (-), the modulator is so biased and driven that its 
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Figure VI-1. Positioning Coils on Tilt Table 



Figure VI-2. Electronic Breadboards for Position Servo 
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Figure VI-3. Block Diagram of Position Servo Electronics 
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Figure VI-4. Position Sensing Characteristic 
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output does not go much below the level it has when the ball is at the quiescent position. This 
ensures the presence of a minimum signal from which to get phase and hence position information 
from that coil. 

C. Dynamic Analysis of Loop 

So far, discussion has been confined to the static properties of the control system. The individual 
elements of the system, shown in Figure VI-3, in general, take a certain time to react, or more 
accurately, the transient behavior is some function of time. Looking at the position sensing first, we 
find that phase angle relates to position. The phase angle is sampled and summed each half cycle so 
that, at 10 kHz, position data lags position by about 25 //seconds. The voltage and current 
amplifiers have no delay networks and, in fact, have several MHz bandwidths. To reduce ripple in 
the position information, ripple filters are employed at the summers and at the modulators. The 
value of their cut-off frequencies are determined by allowable ripple and loop bandwidth. 
Representative values are 300 or 400 Hz for a 100 kHz operating frequency. 

The summers contain, in addition to a ripple filter, a lead network which is required for damping 
the loop. This lead network has associated with it a lag at a frequency of about ten times the lead 
frequency. The exact frequency of the lead network is determined by the gain and bandwidth 
requirements for given load dynamics, the value in this case being about one Hertz. 

The last electronic time constant belongs to the resonated coil which, at 1 0 kHz, has a Q of about 
25, and therefore, a frequency bandwidth of about 400 Hz. The final dynamic consideration is the 
load, in this case, the sphere of mass m. The coil exerts a force f = ki~ for small signals so the change 
in force per change in current is 2 ki v i. We see immediately the square law property of the force, 
which means the overall loop gain is proportional to the current. At a given current, however, the 
force can be represented by k' v i. The transfer function of a mass without viscous force is simply 
X(s) = f(s)/ms^, where X represents distance. This signifies a double lag (-40 dB per decade) passing 
through the point 1/m at 1 radian per second. 

From the above, it is seen that there is no vis ous damping; therefore, it is necessary to provide 
damping in the control loop. This necessity is met by the lead which is placed at the summer. In 
addition to these transient responses of the various elements, each element has a static gain. The 
electronic gains, once set, are fixed, but the position sensing and replusion forces vary. Figures VI-4 
and V-7 show the relations between phase angle and position, and between force, current and 
position. The most variable element is the force which is not only a function of i~ but of the axial 
distance between coil and ball. Therefore, the loop characteristics must be determined for a range of 
static gains. Assuming that the ball stays reasonably well centered in the present loop, this gain 
variation nominally is about 1 0 to 1 , the higher gain being with the ball deliberately decentered 
causing a high coil current to flow. 

A deviation in these dynamic and static gains occurs when the mass is a pendulum bob, as shown in 
Figure VI-5. This deviation is due to a different load transfer function which now is affected by a 
gravity induced restoration force and its attendant natural frequency of vibration. This transfer 
function for the pendulum mass is X(s) = f(s)/m[s“ + g/L] , where L is the pendulum length and g is 
the acceleration due to gravity. 
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Figure VI-5. Force Measuring Apparatus Tilt Table and Pendulum 
The total dc transfer is 

G = 2.86 x .35 x 20 x . 1 x 2. 1 2 = 4.24 volts per volt. 

20 log 1Q 4.24 = + 12.6 dB. 

A Bodie plot of the important frequencies appears in Figure VI-6 for these gains and break 
frequencies. In actual practice, both the resonance frequency and the gain were higher than 
calculated which placed the 0-dB crossover at a somewhat higher frequency and thereby caused an 
underdamped condition. 

D. Position Sensing Considerations 

Experience with the servo showed position sensing as the most critical task. The task here is made 
difficult because (1) the coil excitation changes in amplitude about 3:1 from the servo output, and 
(2) the other coils couple strongly into it. The position sensing logic rejects all but a small amount 
of this disturbance, but the residual approaches the sensitivity level for adequate position sensing. 
Some improvement can be made with the existing circuit; but to incorporate spinup, levitation and 
electromagnetic heating, other approaches offer significant improvement and greater flexibility. For 
instance, position sensing, a small signal function, can be made independent of the force coils, a 
power function, in power, frequency and geometry. Some of the more promising approaches are as 
follows: 
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1 . Synchronous detection of current to eliminate all other frequencies. 

2. Balanced bridge with synchronous detection, single frequency excitation for all 
coils. 

3. Time-sequencing one axis at a time with momentary deactivation of all other fields. 

The position sensing and control device described above has been thoroughly tested and has proven 
itself capable of quickly damping sphere motions and restoring the sphere to the center position 
where it remains at zero relative velocity. The stiffness and damping constant of the loop can be 
easily adjusted. In the experiments carried out under the present contract, initial ball displacements 
were of the order of 1 cm from the center position, corresponding to horizontal accelerations on 
the order of 5 x 10'^g. Thus, control loop stiffness in terms of applied field strengths was 
significantly higher than would be required in a near zero-gravity velocity facility. Testing of these 
devices in a drop tower would be desireable at operating levels more closely approaching those of 
interest to a zero-gravity space facility. Use of minimum restoring forces may be important for 
applications where shape deformation of the suspended mass is to be avoided. 

A five minute color movie was prepared documenting results of position control and damping 
characteristics. 
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E. Spinning a Conductive Sphere 

An aluminum sphere (2.1-cm diameter) and a copper sphere (3.1-cm diameter) were each set into 
rotation while supported by a small ball bearing among the four positioning coils previously made 
for the servo-control demonstration. The four coils were energized at 60 Hz, with the field 
generated by each pair of coils about 90 degrees out of phase with that generated by the other pair. 
The apparatus is shown in Figure Vl-7. With the rotating field at the position of the copper sphere 
approximately 50 oersted, each sphere was easily rotated at a constant rate about one revolution 
per second by balancing the torque generated by the field against that provided by the drag of the 
ball bearing. 

Each sphere was then hung from a 300-cm long wire, which had a torsion constant of 5.8 
dyne-cm/radian, to measure the torque imparted to the sphere by the rotating field. For a field 
strength of about 30 oersteds at a frequency of 220 Hz, the stationary aluminum sphere received a 
torque of 25 dyne-cm and the copper sphere received 1 20 dyne-cm. 

These experiments were later repeated with the ball suspended from a long thin wire. Spin up with 
horizontal forces furnished primarily by the coils was observed to be stable until wire windup was 
encountered. 

A glass sphere (described in Appendix D) was filled with a saturated aqueous solution of sodium 
chloride and also given a torque. With a field strength of 35 oersted at 10 kHz, the torque imparted 
to the assembly was 0.07 dyne-cm. 



Figure VI-7. Solid Copper Sphere Supported on Ball Bearing Centered Among Four Coils which 
Generated a Rotating Magnetic Field 
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VII. GENERAL REQUIREMENTS FOR SPACE EXPERIMENT FACILITIES 


A. Introduction 


Each of the proposals for exploiting the weightless feature of the space environment m metallurgic 
and ceramic processing of course requires an individual study of facility requirements in terms of 
heating power, instrumentation, size of facility, etc. Smce a large number of these suggestions involve 
the handling of a floating molten mass out of contact with crucibles or molds, it is nevertheless 
possible to define, even at this stage, facilities which will be capable of accepting a large number of 
candidate new weightless processing experiments. The following describes the beginning of a study 
to define the physical requirements for a wide range of potential new crucibleless melting and 
solidification experiments in order to define limits to the range of variables which can be handled m 
one or more facilities which may be practicable for early incorporation in the post Apollo program. 

Some of the physical variables which are important m defining the required experimental facilities 
are summarized below 

Process Variable Facility Requirement 


Melting temperature 
Size of batch 
Heating rate 


I Heating power 


Requirements for vacuum or 
controlled atmosphere 

Requirements for pre-meltmg of 
pre-cast specimens or for crucible- 
less mixing and reaction 


Processing chamber and accessories, 
starting material, handling devices 


Cooling rate 

Solidification temperature 

Processing time 

Purity requirements 

Requirements for absence of contact 
with crucible of melt during some 
part of the process 


Position control, free floating volume 
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Requirements for stirring Electromagnetic or electrostatic fields, 

( provision for rotation of melt 

Requirements for shaping or molding 

Likelihood and nature of included Requirements for stirring and bubble detection/ 

bubbles elimination technique 

We must translate some of these rather general process requirements into a definition of the 
facilities required to handle a wide range of processes. The requirements on facilities imposed by 
position control, position sensing, heating, processing time, free floating volume required for the 
specimen, maximum permissible accelerations in terms of shape distortions, etc. are of course all 
related. We first discuss these requirements separately, however, in the course of a detailed 
discussion, the interdependence of these requirements will become obvious and the region of 
mtersection of the various experiment requirements in terms of facilities will, in most cases, define 
the required facilities within fairly specific limits. 

B. Material Properties — General Considerations 

Electromagnetic waves of any frequency (greater than those which might better be called slowly 
varying dc fields) falling upon a conducting material of zero resistivity would behave as light falling 
upon a mirror; all of it would be reflected. If the conducting material has a non-zero resistivity, then 
the electromagnetic waves penetrate the material, loosing energy in the process, and the depth by 
which the amplitude of the wave has dimished to 1/e or 0.37 of its initial value is the attenuation 
distance or the skm depth. This skin depth, 5, depends upon the properties of the material and the 
frequency of the waves as follows* 

1 

6 

where f is the frequency of the wave in cps or Hz, a is the conductivity of the material in 
mhos/meter (one mho/meter is equivalent to a resistivity of 100 ohm-cm), p is the magnetic 
permeability of the material in Newtons/ampere^ (p = 47T x 10'^ Newton/ampere^ if the relative 
permeability is unity), and 8 is the skm depth in meters. The skm depth has been plotted in Figure 
VII-1 for several materials and a wide range of frequencies. Here it may be noted that the skin 
depths achieved m aluminum and copper in work discussed elsewhere m this report ranged from 0.8 
cm at 1 00 Hz to 0.5 mm at 1 0 kHz. 

C. Melting Power Required to Supply Surface Radiated Power Loss 

Table VII-1 presents the melting points of various materials and the total power required for a 
sphere of 10 cm^ surface area to overcome radiation loss at the melting point. This is an important 
factor in selectmg materials for experimentation and in designing apparatus to melt these materials. 
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Figure VII-1. Skm Depth versus Frequency for Several Common Metals and Molten Glass 
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Table VII- L Physical Properties 


Substance 

or 

Element 

T Melt (C) 

r 

T Melt (°K) 

Surface radiation loss for 
10 cm 2 and unity emissivity* 

Ti 

1800 

2073.16 

1.04 x 10 3 

Be 

1283 (Top) 

1556.16 

3.32 x 10 2 

Zr 

1857 

2130.16 

1.17 x 10 3 

W 

3370 

3643.16 

1 00 x 10 3 

Mo 

2620 

2893.16 

3.97 x 10 3 

Cr 

1890 

2163.16 

1.24 x 10 3 

Ta 

2996 

3269.16 

6 47 x 1 0 3 

PbTe 

917 

1190.16 

1.13 x 10 2 

Steel 

1537 (Top) 

1846.16 

6.58 x 10 2 

Ni 

1455 

1728.16 

5.05 x 10 2 

Ga 

29.78 

302.94 

.48 


5.668 x 10' 12 watts/cm 2 (°K) 4 


iffl 











D. Basic Equations for Acceleration and Heating 

The development of a set of parameters to discuss processing and positioning of materials can be 
most easily made by considering the relatively simple system of a single coil and a material 
sphere The same parameters are also the most useful m discussing more complicated systems with 
more than one coil, as will be shown below. 

The processor and designer wishes to relate the power absorbed by the material and the force 
exerted on the material in terms of material properties. For a single loop and a material sphere, the 
force exerted by the loop on the sphere along the axial line of the loop is 


F (dynes) 


_ 3 _ 

50 


2 2 
IT I 


A (y) G (x) 



( 1 ) 


and the eddy current power absorbed by the sphere from the field is 

2 

N (watts) = 3f E p H F 1 (x) 

ct Q 1 


( 2 ) 


where 


1 1 is the current in amperes 

y is the distance along the axis of the loop m terms of loop radii z/Rj, where z is the 
distance along the axial line 

Rj is the loop radius 

R2 is the sphere radius 

x is the ratio of sphere radius to skin depth R2/8 
5 is the skin depth 


A(y) = 


2 4 

(i + y ) 


is the magnetic field intensity from the loop along the axis 



F x (x) 


x (sinh 2 x + sin 2 x) - cosh 2 x + cos 2 x 
cosh 2 x - cos 2 x 


p e - the resistivity in ohm-meters. 


For a single loop 


H = 


4R‘ 


(i+y 2 ) 3 


so that 


N - 


3 7T 
4R 2 


(i + y 2 ) 3 


R 2 ^ *1 


P 1 (X) 


( 3 ) 


(4) 


It can be shown that for more complicated arrangements the equations have the same form so that 
as an example for two equal coaxial loops carrying equal alternating currents 


F = — 


3 

50 


/ ID \ ** 

2 2 I 2 \ 

7T I x G (x) B (y) ^ J 


(5) 


where 


B(y) = 


(i + y 2 ) 3/2 


.. 2 3/2 

(1 y l ) 


2 5/2 - 


(i + y )”' " (1 + y~) 


.2 5/2 
V 


Here the distance z is measured positive upward from the plane of the first loop and zj is measured 
positive downward from the plane of the second loop Then 

_z_ _ \ S 

R 1 R x 7 “ y l 


where S is the distance between the planes of the loops. 


vn-6 



Thus the same parameters used to discuss the single loop and material sphere can be used to discuss 
more complicated systems and the mathematical complexity is reduced to a minimum by 
considering the simplest system 

1. Development of Relations Involving only Matenal and Driving Coil Parameteis 


The parameters y, x, etc., used in the above equations are very suitable for computation purposes, 
but the designer wants to speak of force and power m terms of current, frequency, and the matenal 
parameters. Looking at Equations 1 and 4, it is evident that even for the single loop there are many 
parameters in each equation even with the dimensionless parameters x and y. Thus, the problem of 
speaking coherently about many materials and their behavior under applied fields and plotting 
results is fairly complicated 


The problem was attacked in the following manner. First the distance dependence factor was 
removed. The force curve has a peak close to y = 0.4. Comparisons were made then at A(y = 0 4). 
The loop radius was fixed at Rj =0.75 inches or 1.905 centimeters. This converts Equation 1 to 
F = constant • if G(x) and Equation 4 to N - constant • R 2 P e if F j (x). For positiomng it is 
more convenient to speak about accelerations, g, imparted to the sphere; since the mass of the 
sphere is m = 47rpgR2 where p„ is the density, we have 


S = 


4tt R 


2 "g 


and after putting these into consistent units we have 


( 6 ) 


g (newtons /kilogram) = 


0. 0452 G (x) I 


g 


(7) 


N (watts) = (4.16) . 10 3 p e 1^ R 2 F ;l (x) 


where here p g is the density in kilograms/meters^ 


(B) 


p e is the resistivity in ohm meters, and Ij is the current in amperes. The parameter x is just Rp/5 where 

6 = (V7777 j 

and p is the permeability, o is 1 /p e the conductivity, and v is the frequency of the applied field so that 


x 


= R 9 


(9) 
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The relations of Equations 7, 8, and 9 then form a consistent set for discussing different materials in 
the field of a single loop. 

For a close wound coil we can replace Ij by NIj without an appreciable error, where N is the 
number of turns in the coil. Calculations based on this approximation for a very close wound coil 
have been shown to agree fairly well with experimental results (previous quarterly). 

2. General Behavior of Force and Power with Skin Depth 

It can be seen from Equations 7, 8 and 9 that the following results are obtained. The acceleration 
imparted to a sphere at a fixed distance from the loop depends upon the sphere radius, the square 
of the loop current, the permeability of the sphere material, the conductivity (or resistivity) of the 
sphere material, and the frequency employed in the loop through the complicated function G and 
inversely upon the density of the sphere material. The eddy current power absorbed by the sphere 
from the field depends upon the resistivity, the square of the loop current, the sphere radius and 
also the sphere radius, permeability, conductivity, and the frequency employed through the 
complicated function F j(x). 

For a given material the ratio of power to acceleration N/g is just 

/F l (x) 

constant • i — , - 
\ G (x) 

Figure VII-2, taken from Fromm and Jehn’s article, shows Fj(x), G(x) and Fj(x)/G(x) versus x 



Figure VII-2. Ratio F j (x)/G j (x) and the Functions F j (x) and G(x) 
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Regimes of x and Range of Parameters 


Referring to Figure VII-2, we see that we can speak essentially of three different regimes in x (ratio 
of sphere radius to skin dep th). T hese are the regimes x g l, x from 1 to about 10, and x > 10. For 
a given material, x = R2^7tiiov = constant and frequency is the principal variable, but for 
different materials ranging from dielectrics with p e ~ 10^ to conductors ranging to p c ~10'^ or 
less (ohm meters) with a = 1 /p e , the material parameters influence the parameter x considerably 
Thus for a glass or a semiconductor, a given yf may correspond to very high frequencies, whereas for 
a conductor it will correspond to very low frequencies. Thus, the region x < 1 may range in the 
frequency from tens of cycles/sec to 10^ or more cycles/second depending on the material 


4. 


Small x Power and Force Relations 


It is found that for small x 


(x) « 


4 x 


45+ 2x 


4 x 
45 


G(x) 


8 x 

315 + 14 x^ 


8 x 
315 


(very small x) 


(very small x) 


( 10 ) 


(ID 


and that at very small x, Fj/G « 3.5. This is shown in Figure VII-2 where below x ~ 2.5 the curve 
of F|/G flattens out at the value 3.5. 

For small x (actually it is found that this relation is fairly good out to x w 2.5), it is found that the 
power to acceleration ratio is 

- = 3.22 Rp p (12) 

g 2 r g 'e 

with N/g in watts/(newtons/kgm) 

Pgingm/cm 3 ] for 

convenience 

p e m microhm-cm J 

For a specific material under a fixed force we can see that N = constant * R 2 so that, for a given 
acceleration and a specific material, the absorbed eddy current power depends only on the sphere 
radius. For each material, then, the absorbed power will depend only on the sphere radius at small 
x, and the plot of N vs R? on log-log paper will be a straight line. 

5. Positioning Control at Small x 

Figure VII-3 shows the eddy current power absorbed by a free floating specimen when it is 
accelerated by an eddy current positioning device at a rate 10‘^g (10 3 cm sec ). This acceleration 
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is characteristic of the lower limit which may be achievable in lower altitude earth orbitmg facilities 
unless compensation for decelleration due to air drag and gravity gradients are provided. The power 
requirement depends upon the sphere mass, density, electrical conductivity, and driving frequency. 
It also depends upon the magnetic permeability of the specimen, but our main interest is in molten 
materials which will be above the Curie temperature so that we may assume unit permeability. We 
see that the power is linearly proportional to the sphere radius and depends drastically upon the 
conductivity of the material. For example, minimum position control power for a one centimeter 
aluminum or molten glass sphere at an acceleration of 1 0"^g ranges from 1 CT~* watts to 1 00 watts. 
For molten glass of relatively high conductivity, the corresponding power required for position 
control drops to 10 milliwatts. The corresponding frequencies required for position control range 
from the low audio range for metals to tens of megahertz for the molten glass of poorest 
conductivity. 


For relatively small experiment vehicle carriers, such as the proposed dry work shop, astronaut 
body motions may occasionally impart accelerations to the vehicle approaching ICT^g for short time 
periods. For most processes, sufficient free volume can be provided surrounding the specimen being 
processed so that it can remain free floating without seeing accelerations of this magnitude. Since 
these accelerations will average out to zero over a time period on the order of a minute or less, it is 
highly likely that rigid position control with respect to the vehicle will not be required. At least for 
short processmg times, the operation of such accelerating forces could easily be inhibited until 
process completion. If it is desired to provide rigid positioning of free floating object with respect to 
the vehicle during accelerations as high as lfH'g, we can refer to Figure VII-4 for the required 
positioning powers. For the metals in the range up to ten centimeter radius, the positioning powers 
are essentially negligible (less than one watt). For high conductivity molten glass, the power 
requirement ranges up to only ten watts. For poorly conducting molten glass, required positioning 
powers can reach the multikilowatt region and hence it appears that eddy current position control is 
not a likely candidate m this extreme case unless simultaneous RF heating is used 


Referring to Equations 7, 8 and 9, and using the small x approximation, we see that for x < 1 or 
R 2 /5 < 1 



constant • 


4 

x 


P 


g 


N 


4.16 


103 Pe h 



constant • 



R, 


4 

x 


and for a fixed material, x = constant R 2 so that 




N oc 


T 2 n 5 

h P e R 2 


2 

V 


for small x. Small x (x < 1) means < 1 so that the skin depth of penetration exceeds the 
sphere radius R^. 

For the same material, then, for small x 
2 

g oc v 

Kocj / 2 

Thus, for material conductivities, sphere radii, and driving frequencies such that the skin depth of 
penetration exceeds the sphere radius, a change in driving frequency will cause the acceleration and 
power dissipation to vary in the same proportion. Increasing the frequency so as to double the 
acceleration will also double the power dissipation. 


6. Positioning at Large x 

Curves showing the manner in which the acceleration and power dissipation vary with frequency are 
shown in Figure VII-2. These curves show that for frequencies and conductivities such that the skin 
depth is a small fraction of the sphere radius, the power dissipation increases with frequency more 
rapidly than acceleration. For the metals, positionmg with negligible heating at low frequency is 
relatively simple. For a very poor conductor such as molten glass, an increase in frequency to at 
least the high audio range and even into the tens of megahertz range is required, hence, eddy current 
heating cannot be avoided. We have displayed the dependence of power dissipation, normalized per 
unit sphere surface area, as a function of driving frequency for metals and semiconductors in Figure 
VII-5. Figure VII-6 shows the range of frequencies of interest for positioning of metals and 
semiconductors. 


For the molten metals of Figure VII-5, we can no longer assume, for practical frequencies, the 
regime of small x. Calculations for this plot are based on the full expression for F j(x) and G(x), 
case x ^ 1 given earlier. 


E. Frequency Choice for Facilities 

Here we see that the choice of frequency for a positioning facility can be made conveniently in the 
audio region if the specimen is heated by some means other than eddy currents. In practice, we can 
select a frequency somewhat arbitrarily, such as ten kilohertz, for winch requirements for electronic 
circuit design are extremely simple. The actual positioning forces obtained at a given driving field 
strength will then of course depend upon the resistivity and size of the material sphere being 
processed In the servo positionmg device which we are developing, the exact restoring force per 
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Figure VII-6. Acceleration as a Function of Dnving Frequency for Metals and Semiconductors 




unit position error is not particularly critical. Moreover, the servo loop gam can be varied quite 
easily, if desired, for different specimens by turning a knob to a precalculated position. In practice, 
such an adjustment would probably not be made unless the material resistivity varied extremely 
from one specimen to the next. 

For molten glass (Figure VII-7), we see that the frequency required to obtain positioning 
accelerations of 10'^ g range from one kilohertz at 1,000 ampere turns and a resistivity 10"^ ohm 
meters to 100 megahertz for the case of a one ampere turn driving coil and glass of resistivity one 
ohm meter. It thus appears that, although eddy current positioning is feasible for molten glass, it 
may not be the optimum method in practice, especially for glass of poor conductivity. We have 
previously showed the skin depth as a function of frequency in Figure VII-I. Depending upon the 
frequency and conductivity, we see that the skin depth can range over 7 decades of variation, thus 
the ratio of skin depth to sphere radius can vary from values much less than unity to values much 
greater than unity for cases which will be of interest. This is in contrast to work in terrestrial 
levitation experiments where skm depths greater than the sphere radius are of little interest since 
they, in general, do not provide sufficient positioning force to levitate the material in the one g 
environment. The greater freedom of choice of frequency and skin depth for the case of matenals 
which are processed in the zero gravity environment allows for eddy current heating which can be 
adjusted independently from the requirement for positioning. 

F. Heating and Temperature Control 

The minimum hearing power which must be provided is that necessary to furnish surface radiation 
loss from the specimen If a controlled atmosphere is utilized, additional power will be conducted 
away by the gas. In this connection it must be noted, however, that normal convection due to 
buoyancy of the heated gas will be absent and that, under these conditions, the gas will act as a 
good insulator. For many matenals, including most metals, the choice of hearing methods can be 
made independently of the choice for position control. For matenals of very high resistivity such as 
glasses, the electromagnetic positioning can give nse to eddy current heatmg dissipations which can 
become significant so that it may be advantageous to consider RF heating, at least for things such as 
molten glass having reasonably good conductivity. One might be led to think, conversely, that the 
choice of eddy current or RF heating will lead to unavoidable positioning forces. However, if RF 
heating is provided by pairs of coils symmetrically disposed on either side of the specimen, no net 
translational force will be applied to the specimen when the coils are equally excited. With such an 
arrangement, if it is desired to impart a position control force, small differential excitations can be 
introduced mto the coil excitations through a servo loop actuated by position or velocity errors in 
such a way that the total power dissipation in the specimen is unchanged. The specimen 
temperature can be controlled by raising or lowering the total coil excitation power. 

Figure VII-8 simply illustrates the surface radiation loss, or minimum hearing power required to 
achieve the meltmg or transition temperature for several materials of interest. By way of illustrating 
the relation between the heating power requirement, if this is provided by eddy currents, and the 
eddy current positioning requirements discussed earlier, we may consider the example of PbTe 
shown in Figure VII-8 This requires a minimum power dissipation in a one centimeter radius sphere 
of 130 watts. Referring back to Figures VII-5 and VII-6, we see that this requires a driving 
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Figure VII-7. Power as a Function of Dnving Frequency for Molten Glass 





Figure VII-8. Minimum Heating Power Required to Achieve Meltmg Temperature for Several 
Metals of Interest 

frequency of 600 kilohertz at 300 ampere turns, corresponding to a maximum acceleration of 0.04 
g if all of the power is furnished by a single coil. Any lower value of position control acceleration 
can be provided by sharing the driving power between two coils, one on either side of the specimen 

G. Position Sensing 

A number of position sensing schemes are available which can be used for position control or 
handling of free floating solid or liquid materials. Besides the eddy current position sensor which we 
have developed and applied to sensmg of metallic objects, one may also consider bolometers or 
other heat sensmg detectors or changes in electrostatic capacitance. The latter would be a natural 
solution m the case of an electrostatic position control device. 

Another possibility is optical sensing, either by emitted radiation from a hot object by interception 
of optical light beams, or visual or photographic observations It is expected that visual observations 
of specimens during processing should be provided wherever possible, particularly m pilot 
experiments. 

Any of the above position sensmg methods, except photographic, can be used to actuate a position 
control servo m the same manner as is done m our demonstration electromagnetic device. It is quite 
likely that none of the above methods will be applicable to all processes, but each will have its own 
domain. It appears that these domains overlap. For many processes, either the eddy current position 
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sensor or the bolometer position sensor will be adequate. For samples of very high resistivity which 
must be positioned at temperatures too low for detection of emitted radiation, the electrostatic or 
optical methods can be considered. Normally, however, positioning of these specimens can perhaps 
be performed mechanically when they are removed from the processing chamber after cooling. 
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APPENDIX A 


MEASUREMENT OF SMALL FORCES 

The forces of repulsion measured in 1969 were initially of the order of hundreds and thousands of 
dynes. After a series of refinements were made, forces as low as five dynes were measured acting 
upon an aluminim sphere of 13.5 gram mass (force/weight 
intention of attaining force weight ratios closer to 10"^ as our 
recently reviewed more closely. The techniques examined were 

1 . Sphere free to roll on a flat surface. 

2. Simple pendulum. 

3 Sphere suspended in a hquid. 

4. Torsion balance. 

5. Electromagnetic levitation. 

In each case it was assumed that the mass being pushed by the divergent electromagnetic fields was 
a 30-gram steel sphere with a 1-cm radius. 

1. Sphere Free to Roll on a Flat Surface 

It had been planned to employ a plate glass tilt table as soon as possible m the measurement of 
small forces, but it was soon evident that, for speeds and forces with which we wanted to work 
(e.g., 0.1 cm/sec and lO^g), the presence of dust presented a problem. Also, it turned out to be 
desirable to use a good conductor such as aluminum to facilitate position control at low frequencies 
and a sufficiently high quality aluminum sphere was not readily available. The use of the glass plate 
was set aside, then, and we turned to the simple pendulum described in the quarterly report of 
November 15, 1969. This rolling technique has now been reconsidered m more detail and the 
following limitations have been recognized. 

a. Dust - In a laboratory in winch no special attempt is made to clean or filter the air, a 
rough estimate of the force required to start a perfect sphere rolling on a perfect 
surface might be as high as 1 00 dynes. Use of a moderately clean-room environment 
could reduce it to something more like 0.1 to 10 dynes The dust particle density 
requirement for several classes of clean rooms is given'in Figure A-l . Sampling of the 
air in such a room must produce data points which lie below the proper line in this 
figure. Assume that the surface upon which the sphere is to roll is in a closed 
container (such as our Plexiglas pendulum box) which has been cleaned such that it 
meets the requirements of a class 100,000 clean room. Assume, also, that all the 
dust in the air m this container up to one inch above the surface, settles upon the 
surface. The remainder then clings to the walls or remains suspended m the air 
during the course of the experiment. The distribution of particles corresponding to 
the line for a class 100,000 clean room is 


ratio equaled 4x1 0"^). With the 
goal, the techniques available were 
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where n is the particle density (ft“°) and s is the particle size (microns). The particle 
density on the surface is then 


N 


■W) 


1A 4. 3 -2.3 

10 s 


where N is expressed in ffi^. Let us now estimate to what size particle we must go, 
starting from the largest, to find a surface density sufficiently large, on the average, 
to prevent a sphere from rolling freely. 


A sphere can be prevented from rolling freely if it is surrounded by three particles of 
equal size as shown in Figure A-2. The surface density corresponding to such a 
condition is A ' 1 (ft" z ), where A is defined in the figure. Thus, setting N = A' 1 and 
arranging for R to be in cm and s in microns: 


s = [ 2. 4 x lcf 4 R] ~ 1 ' 3 


If R = 1cm: 


s = 0. 002 micron. 

The smallest force required to push the sphere away from the three trapping 
particles is approximately (for s < < R): 


F = 


W 

2 



and for a 30 gm sphere 


F = 6 dynes 

Because the class 100,000 line m Figure A-l is an upper limit for any particle count 
m the clean area and because of the inexact nature of this calculation, the starting 
force might better be given as probably in the range between 0.1 and 10 dynes 
depending upon the dust distribution on the surface. 

b. Sphencity - The best steel balls have a tolerance on then sphericity of ± 0 000010 
inch. If an out-of-round condition of 0.000010 inch is assumed to be in the form of 
a flat spot, a force of 200 dynes will be needed to initiate rolling In actuality, the 
out-of-roundness is not found to be in such an extreme and simple condition and a 
more realistic mitiatmg force is probably in the range of 20 to 2 dynes. As in the 
case of dust, any out-of-round condition would affect the path of a slowly rolling 
sphere. 
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Figure A-l. Maximum Density of all Particles with Diameter Greater than s for Several Classes of 
Clean Room 
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Figure A-2. Preventing a Sphere from Rolling Freely 
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c. Fnction - Data for steel spheres rolling on a glass surface have not been located but 
from data on similar materials it seems that a coefficient of rolling fnction of 0.001 
may be expected with no possibility of improvement beyond 0.0005 This provides 
a force/weight ratio of 10“^ to 0.5 x 10’^ and a force of 30 to 15 dynes for a 30-gram 
mass. 

2 Simple Pendulum 

A simple pendulum one meter long with its position observed through a microscope or telescope, as 
done previously, or via an optical lever (reflected light beam) may be just barely capable of 
measuring force/weight ratios of 10"'\ The feasibility depends upon the quality of the optical 
components used and largely upon the degree to which the effects of vibration can be minimized. 
Standard microscopy objectives cannot be used in the vicinity of the suspended sphere if conductive 
materials are to be prevented from distorting the electromagnetic field. If this technique were to be 
pursued, probably the layout of the apparatus would be as shown in Figure A-3. 

3. Sphere in Liquid 

To give a mass of material complete freedom to move in two dimensions it could be suspended by a 
float in a non-conductmg liquid shown m Figure A-4. The equation of motion for small v is 


m — + 3ir 6 vp = F 



Figure A-3. Layout of Apparatus for Measurement of Small Forces with a Minimum of 
Electromagnetic Field Distortion 
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where v is the velocity of sphere relative to the fluid, 5 is the diameter of the sphere, p is the 
dynamic viscosity of the fluid, and F is the external force applied to the sphere. The second term of 
this equation is the drag force as given by Stake’s law which holds for Reynolds’ numbers less than 
0.5 and 

Reynold's number = 

where p is the density of the fluid. One approach to this technique would be to observe the motion 
of the mass immediately after the initiation of the electromagnetic field, while v is very small and 
the viscous drag force exerted by the liquid upon the sphere is small compared to the applied force. 
During this period, the applied electromagnetic force would be accelerating the mass as if it were 
floating freely m space. However, the several requirements tended to be mutually contradictory. 

For example, let it be required that a constant force F be applied to the initially stationary sphere 
and, at the end of the period of observation during which a distance s has been traversed, the drag 
force is still 1 0 times less than F. That is, 

3 ir 6 vjiS 0. IF 

If F = m x 10'~* and g = 0.3 dyne for m = 30 gm in water 


v < 0. 003 cm/sec 




for which the Reynold’s number is § 0.006. The distance s is 0.0005 cm which would be traversed 
during a period of 0.3 second. Measuring these last two quantities with reasonable accuracy is 
possible, but not without considerable difficulty and sophisticated specialized apparatus. 

Observations can be made readily when the applied force equals the drag force. If the terminal 
speed thus attained produces a Reynold’s number less than 0.5, then Stokes’ Law would permit 
direct calculation of the average applied force during the period of observation. The drag force 
exerted upon a 2-cm diameter sphere in water, traveling at 0.05 cm/sec (for which Reynold’s 
number is 0.01), is 0.1 dyne. The time to travel each of several 1-mm observation lengths of its path 
would be 20 seconds. The flotation device would add some drag which would have to be accounted 
for. 

4 Torsion Balance 

A highly sensitive torsion balance technique has long been used m our laboratory for measuring the 
force molecular beams exert upon a surface, and this technique should be adaptable to the needs of 
this program. For example, a tungsten filament 0.05 mm in diameter, 10 cm long, will provide a 
torsion constant of 0.13 dyne-cm/radian. If the mass upon which the electromagnetic field is acting 
is suspended on a 7.5-cm radius from the bottom of the tungsten filament, a force of 0.06 dyne will 
move the mass horizontally 2.5 mm (i.e., 0.02 dyne/mm). Such a force acting upon our 13.5 gm 
aluminum mass, for example, would produce an acceleration of 0.5 x 10'^g and upon a 30-gram 
mass, an acceleration of 0.2 x 10'^g. The basis of this calculation and an estimation of the 
uncertainties involved are given below. The torque provided by a torsion balance is given by 


2 L 

where 6 is the angle of rotation, E is the shear modulus of the filament, R is the radius of the 
filament, and L is the length of the filament If the mass is at the end of an arm of length X from the 
filament’s axis, then the force acting upon the mass is 


7T 9 E R 4 

F - 

2 LA 

This equation is adequate for designing a torsion balance to measure small forces, but it does not 
lend itself to' accurate force measurement because of the term. If the filament’s diameter and its 
maximum uncertainty were 0.001 + 0.00005 inch, the uncertainty in R would be 5 percent and the 
uncertainty in R^ would be 20 percent. Add to this a total uncertainty of 5 to 1 0 percent for all the 
other terms and the force would be known within + 25 to 30 percent The calibration of this 
torsion technique is much improved if the period of a torsional oscillation of the system were 
measured with and without a mass of known moment of inertia. A cylindrical rod is very useful for 
this purpose and the force is given by 

2 

_ it m 6 i 
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where m is the mass of the rod, A is the length of the rod, t r is the period of oscillation of the 
balance with the rod lying across it, and t^ is the period of the balance alone. Cumulative errors in 
this case should not exceed 10 percent, especially if m and A are large. 

The tensile strength of a 0.05-mm diameter tungsten filament has been estimated and, according to 
this estimation, it should be capable of suspending a 900-gram mass. For the purpose of this small 
force measurement, no more than about 20 percent of this limit would be required to suspend and 
balance a 30-gram mass, and an additional 10 to 20 percent would be used for the extra calibration 
rod. 

A variation of this torsion balance which might be useful in demonstrating the dynamic behavior of 
a mass under the influence of small forces, much in the same manner that a long simple pendulum 
was used in our filmed demonstration with larger forces, is a mobile made up of two torsion 
balances, one suspended from the other. If sensitivity constants of the order of 0.002 dyne/mm are 
obtained by proper choice of parameters, the test mass would appear virtually free over a 1 cm x 1 
cm range. A disadvantage to this mobile technique lies in the fact that the dynamics of the sphere 
would be tied to the moments of inertia of the system and not solely to its own mass. 

5 Electromagnetic Levitation 

Electromagnetic levitation usually is accomplished m such a manner that a potential well is formed 
at the top of the levitating field to keep molten levitated material centered in the field. We believe it 
is possible to produce a long level potential trough through use of a rectangularly shaped coil such 
that material suspended by it would behave as if it were at the end of an extremely long pendulum 
or on a slippery surface with an extremely low coefficient of friction. The degree to which this 
might be useful in small force measurement has not been investigated. Interactions between the 
levitating force field and the field creating the force to be measured may interfere with the 
measurement. 
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APPENDIX B 


DIELECTRIC POSITIONING 


The theoretical work on dielectric spheres in electromagnetic fields has not been advanced to the 
level where we have studied conducting spheres. An order of magnitude calculation has been made 
for the magnitude of the electric fields necessary to exert a few dynes of force on a Pyrex glass 
sphere of radius a = 0.4125 in. Consider a dielectric sphere in a time harmonic electric field. To the 
first order an induced polarization 

~ 1 

n e + 2 


P = 4ir e a 
o 


E 


A 

k 


A 

will be induced in the sphere by a field E cos co t k, where 
P is the induced polanzation 
e Q is the permittivity of free space 
P e is the relative permittivity of the dielectric. 

To the first order the force acting on the dielectric sphere is F = (P • V ) E, and this is 


F 


4a e 

o 



~ 1 
7? e + 2 


E ^ k 
c> z 


We shall take a field whose gradient is given by dE/dz = E/3 (volts/cm/meter) or E/3 x 10^ 
volts/meter^, where we have expressed E in volts/cm for convenience, since we are using apparatus 
of this size, so that 


F = 4ir e 
z o 



R' 


3 E 


x 10 


Usmg 


R _2 = 0 4125 in. = 0 0104775 meters 

^ = 51 for Pyrex glass 

e Q = 8.854 x 10'^ farads/meter 
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we have 


2 -10 . 

F = 2.46 E x 10 dynes. 

Therefore, to get a force of 2.46 dynes with a field having the above gradient, we would need 
E = 10~* volts/meter or 10 3 volts/cm. 

This result shows that dielectric positioning of a small Pyrex glass sphere requires fields of the order 
of 10 3 volts/cm to get a few dynes of force. For positioning in the weightless state, this is not an 
unreasonable field to generate. 
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APPENDIX C 


SUPERPOSITION OF FORCES AND DISTRIBUTIONS CONTRIBUTED BY TWO COILS 

Further insight was obtained into the addition of repulsion forces when more than a single control 
coil is excited. The body force due to eddy currents in the conducting sphere is proportional to the 
magnetic field strengths within the sphere, which determine the amplitude of the eddy currents, and 
the magnetic field gradient. The eddy current strength j j due to a single coil is proportional to the 
current ij flowing in that coil times a geometrical factor which can be written as f(yj); yj is the 
coil-to-sphere distance. Thus, j ] = ij f(yj). For two coils, j = ij f(y j) + i 2 f(y 2 )- The gradient of the 
field is proportional to the current in a given coil times another geometrical factor which may be 
denoted as g(y). Thus, 


dB 

dy 



g (y x > + i 2 


g (y 2 ) 


for two coils. When two opposing coils are excited, the total field at the sphere position depends 
not only on the amplitudes of the currents flowing in the two coils, but also on their relative phase. 
If we consider that the driving frequencies of the opposing coils are the same (such a servo design is 
possible, but has not yet been implemented in the laboratory), then the total field at the sphere 
position will depend on the constant phase differences between the two coil currents as well as the 
individual current amplitudes. 


If we write, as an abbreviation, g(yj) = Gj cos cot and ij f(y j) = Fj cos cot, etc., we can write the 
total eddy current force acting on the floating mass as j grad B = [Fj cos cot + F 2 cos (cot + $)] 
[Gj cos cot - G 2 cos (cot + d>] , where we have assumed that the current in driving coil no. 2 leads by 
a phase $ the current in driving coil no. 1 . If we take the time average of this expression, we obtain 
the average body force: 


1 COS ^ 

<F> -I (F 1 G 1- F 2 G 2) + — < F 2 G 1- F 1 G 2> 

If we consider the total power dissipation m the floating mass, we can write: 
P ~ j 2 = [ F cos w t + F 2 cos ( u t + $) ] 2 


for the instantaneous power and the average power becomes 


<P>- 


2 2 
F, + Fn 
1 2 


2 


+ F F cos $ 
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We see, as the phase difference of the current in the two coils varies between zero and 1 80°, the 
average power varies between the limits: 


2 P 

<V F 2> 


< F 1 + F 2 > 


The average force over this same range of phase difference varies between the limits* 


< F 1 + F 2 > < G 1 


G 2 > r < F 1 - F 2 > < G 1 + G 2 » 

< F ^ 


If we plot the average power as a function of phase difference we obtain the curve shown in Figure 
C-l. If we further make the assumption that the sphere is centrally located between two coils of 
symmetrical construction so that the difference between F j and F 2 and between Gj and G 2 is due 
to different current amplitudes, we can write 



The average force as a function of phase difference gives the plot shown in Figure C-l. 


It is thus apparent that considerable flexibility exists with this scheme for varying force and power 
dissipation individually. For example, if the sphere is centrally located between two coils with equal 
excitation and it is desired to apply a body force to the mass without changing the power 
dissipation in the mass, we can write 


AP = i 1 Ai l + i 2 Ai 2 + cos $ (i 1 Ai 2 + i g Ai^) 



Figure C-l . Average Power as a Function of Phase Difference 



where Ai is a current increment. 


For 

*1 = ig’ AP = i 1 + Ai 2 ) ( 1 + cos f>) 

If we select A ig = -Aij, we see that there is no change in power dissipation, but a body force is 
obtained. 
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APPENDIX D 


ELECTROSTATIC POSITIONING OF A DIELECTRIC 

A thin walled, hollow glass sphere 3 cm in diameter, hanging on a fine wire about 300 cm long, was 
filled with tap water and placed among the electrodes shown in Figure D-l. The flat plates were 
four cm apart, the tip of the pointed electrode was 3 cm from each plate, and the sphere surface 
was 2 mm from the tip. (The dielectric constant of water and glass at room temperature are 80 and 
5, respectively.) When 3000 volts were applied between the pointed electrode and the two flat 
plates, the uncharged sphere was repelled from the tip by approximately 2 mm which mdicates a 
force of 10 dynes was acting upon it. Increasing the voltage to 5000 volts increased this deflection 
only slightly. One possible problem which can arise in the use of electrostatics was graphically 
shown when the sphere was made to swing such that it touched or very nearly touched the tip of 
the pointed electrode. It acquired a charge with each swing which by simple electrostatic repulsion 
caused the extent of each swing to become greater until finally, after 5 to 8 successively larger 
swings, the sphere was strongly attracted to one of the two plates, where it remained. 

Pointed probes with spacings similar to those used above were used m place of the plates as shown 
in Figure D-2 without any obvious deflection for voltages up to 5000 volts. 



Figure D-l. Hollow Glass Sphere Filled with Water for Test of Electrostatic Repulsion with One 
Pointed Electrode 
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Figure D-2. 


Hollow Glass Sphere Filled with Water for Test of Electrostatic Repulsion with Four 
Pointed Electrodes 
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APPENDIX E 


REFLECTED IMPEDANCE DUE TO EDDY CURRENTS 


1. Introduction 

In this appendix we shall develop more quantitatively the equations governing measurement of the 
position of a floating mass by electrical means. The simplest method, which has been described m 
this report and which served as a basis for the servo loop demonstration experiments, is to 
continuously monitor the fields generated by induced eddy currents in the floating mass. These 
fields will generate emfis in the control loops used to position the mass, and these emfis can be 
detected as equivalent changes in resistance and inductance of the coils. The induction effect 
follows from the fact that the eddy currents have, in general, a different phase than the coil 
excitation currents. In the mam body of the report, circuits were described which detect these 
impedance changes independently of the level of coil excitation. In the following we develop the 
equations for the change in driving coil impedance due to the eddy currents. 

Consider a single loop and a conducting sphere and let O at the sphere center be the origin of 
coordinates. Looking down on the subsystem, we have the illustration 



We may write the total vector potential A Q outside the sphere as 


A 

o 



+ A 


e 


where A c is the vector potential due to the exciting coil and A e is the vector potential due to eddy 
currents. 

If we make the Fromm and Jehn approximation that “a current in a circular loop produces 
approximately the same field as a magnetic field uniform in space and sinusoidally oscillating with 
time” then 
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where 
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and I + 1/ 2 (v) denotes a Bessel function of order ~ Yz. 

The flux through the loop contributed by the eddy currents is just 

$ = 



where S is a surface capping the loop. Using Stoke’s Law, 
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where i is the path of the loop. 


Then 
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R = R^ cos 0 t + R^ sin 0 § 


and 


dR = dl = R. 
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where the time dependence eJ w *has been added and 
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is the induced voltage in the loop due to eddy currents in phasor form. Now letting B 0 be the field 
due to the loop at the sphere center, if the sphere were not there 


B = 
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and 
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is the induced voltage in the exciting loop due to the flux arising from eddy currents. 
Now D = Re (D) +3 m (D) = U + j W 


E = - 
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which is the induced electromotive force rise in the circuit. The voltage drop is -E. 
so we have 
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as the reflected resistance and inductance, respectively. 


In the high frequency limit as co -*• «> 
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Thus we have a decrease m the effective inductance and an increase in the effective resistance 

For a loop of N terms E = -N d3>/dt and if we assume that the loop is close wound so that the field 
B 0 is approximately 
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then 
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2 Conclusion 

Limited expenmental tests were made for the effective change in inductance as the position of a 
conducting aluminum sphere was varied. Although the results varied with y m the way indicated 
from the formula above, a discrepancy of approximately a factor 2 was observed This is considered 
to lie well outside of expenmental errors and is thought to be attributable to the rather strong field 
gradient which was not taken into account in computing the dipole moment due to eddy currents. 
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APPENDIX F 


ELECTRIC DISPLACEMENT CURRENT IN DIELECTRICS 

When considering the properties of conductors and dielectrics, we must consider the form of 
Maxwell’s equations within a material medium. The equation relating magnetic field intensity and 
current for linear, isotropic, nonpermeable media can be written as 


V x (n H) * u J + p k € 3 E/at 

U U U fc) U 


or 


VxB~/x J + fi k e ^E/^t 
o o e o 

\ \ 

Conduction Displacement 

Current Current 

in the M.K.S. system of units, where /i = 47r x 10*^ henry/meter and is the permeability of free 
space, e Q = 8.854 x 10' 1 ^ fared/meter and is the permittivity of free space, and tc Q is the relative 
permittivity of the linear, isotropic, homogenous, nonpermeable medium. 

In a good conductor 


<‘o fl> % lc e £ o 
and in a poor conductor 

"c 7 * Wo ai/3t 


so that ratio 


Q = 


sE/at 

a J 
0 


is a quantity which expresses the importance of conduction currents or displacement currents in a 
material. For a time harmonic field E = E e JC0t , and 9 E /St = j co E, so with J = aE, the ratio Q is 
just K e e 0 co/a. 
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The quantity 1/Q is the ratio of conduction to displacement currents. In a good conductor like 
copper, 1/Q will Tange from lO^tolO^as the frequency increases from 10 to 10^ cps, while in a 
cold glass, 1/Q will range from 10“^ to 10'^ over the same range of frequencies. Figure F-l shows 
the range of l/Q’s possible in the frequency range 10^ to 10^ cps for a molten glass. Here the range 
of conductivities is 1 mho/meter to 10^ mho/meter, and taking a typical permittivity of K e = 5 
(data show k & can vary from 3 to 18), the range of cr//c e e 0 is 2.2 6x 10^ to 2.26 x 10^ and values 
of 1/Q within this rectangular region will range from about 1.59 x 10^ to 1.59 x 10^. it can be 
seen that as we go up in frequency, the quantity 1/Q decreases, since here displacement currents are 
becoming more and more important. We can also conclude that within this frequency band 
conduction currents still dominate appreciably. For intermediate temperatures, the conduction and 
displacement currents in glass will become of comparable magnitude. 





Figure F-l. Plot of p/« e e 0 = (27r/Q) f with Dashed Lmes a Molten Glass Region 
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APPENDIX G 


CO-VARIANT SOLUTION OF ELECTROMAGNETIC EDDY CURRENT PROBLEM 

In the previous quarterly report, expressions were derived for the magnetic field and current density 
inside a conducting sphere immersed within an otherwise uniform oscillating magnetic field. 
Because of the complexity of the mathematical expressions, a computer was used to calculate the 
distribution of fields and forces. Even these solutions, however, are highly idealized m that practical 
applications require utilization of nonuniform fields in order to obtain position sensing and control 
forces. An essential part of the complexity of the problem is related to the requirement for the use 
of spherical polar coordinates and the resulting complication m such mathematical expressions as 
the Laplacian V ^ and the curl. Because of the fact that expressions such as these have the same 
form independent of coordinate system when a covariant approach is taken to the solutions of the 
electromagnetic field and forces problem, this approach offers a hope for simplification of more 
complex problems such as those of nonuniform fields. 

We will denote the components of the magnetic vector potential as A^, m = 1,2,3. The 
electromagnetic field components are represented by the antisymmetric tensor 0 ^, where the 
indices m and v range over the values 1,2,3. The identification is made -that 023 ~ ^13 ~ %> 

^1 2 = -B r Other quantities which we will require in the succeeding discussion are the following 

The “metric tensor” components g^ v and their inverses g^ v relate the expression for arc 
length to the coordinate differentials. In spherical polar coordinates, these quantities are g^ = 1 , 
g 22 = r ^’ §33 = r 2 sin 2 0. The non-diagonal terms vanish 

The determinant of the metric tensor, g, is in this case 
g = r sm 0 

The coordinates r, 8, <E>, are denoted as xj, X2, X3. 

The electrical conductivity is denoted by cr and the magnetic permeability by ix. 

Maxwell’s equation for the curl of the magnetic field is then written simply as 

where is the contravanant field tensor density^ <fP v , and is the'contravariant current density. 

If we assume that the only non vamshing component of A is A 3 and that this depends only on the 
coordinates xj and X 2 , the above expression reduces to 
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Introducing the more usual notation r , 0, assuming that A depends on the time as exp (icot), and 

introducing the expression for skin depth pace = 1/6, we obtain 
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If we introduce % = -v/i r/5, this becomes 
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If we assume that A can be written as the product of separate functions of r and a function of0,A^ 
= R(r) © (0), we obtain the pair of equations 
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where a. is the separation constant. The separation constant satisfying the boundary condition of a 
uniform applied field gives a= -2. This allows the angular equation to be integrated as 

2 

0 ( 0 ) = const x sm 0 
The solution of the radiation equation is 
R (r) = 4 Jy 2 {£ ) 

where J is a Bessel function of the first kind. The complete solution inside the sphere is then 
A g = C sin 2 0 ^ J g/2 (g ) 
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The corresponding unitary magnetic vector potential is 
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in accordance with the solution found by Smythe. 

The expressions for the magnetic field components and the electromagnetic force can be written in 
very simple form in terms of the covariant component A3 as follows: 
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The simplicity and symmetry of these expressions is evident and offers hope for simplification over 
the expressions in unitary form for more complicated problems, such as that of the nonuniform 
field. It must be remembered that the more familiar unitary force component per unit volume is 
related to the covanant force component by the expression 
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